On Quantum Bosonic Solids and Bosonic Superfluids 
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We review the nature of superfluid ground states and the universality of their properties with 
emphasis to Bose Einstein Condensate systems in atomic physics. We then study the superfluid 
Mott transition in such systems. We find that there could be two types of Mott transitions and 
phases. One of them was described long ago and corresponds to suppression of Josephson tunneling 
within superfluids sitting at each well. On the other hand, the conditions of optical lattice BEC 
experiments are such that either the coherence length is longer than the interwell separation, or 
there is too small a number of bosons per well. This vitiates the existence of a superfluid order 
parameter within a well, and therefore of Josephson tunneling between wells. Under such conditions, 
there is a transition to a Mott phase which corresponds to suppression of individual boson tunneling 
among wells. This last transition is in general discontinuous and can happen for incommensurate 
values of bosons per site. If the coherence length is small enough and the number of bosons per site 
large enough, the transition studied in the earlier work will happen. 



I. INTRODUCTION 

We have recently found Ij that bosonic systems in 
optical lattices can have discontinuous transitions from a 
superfluid phase to a Mott insulating phase irrespective 
of conimensurability of the number of bosons and the 
number of lattice sites. This is contrary to expectations 
from the early theoretical work 0. The discrepancy is 
not due to incorrectness in the early work, but due to 
that, depending on the physical conditions, the nature of 
the transition and the insulating phase can be different. 
Specifically for a small number of bosons per site and/or 
for superfluids with coherence length larger than the well 
size, the transition studied in the early work is impossible 
, as it presupposes a superfluid order parameter within 
each well. 

The nature of the Mott phase of the early work is 
that of wells of superfluid, where the large repulsion pre- 
vents Josephson tunneling much as it happens in Joseph- 
son junctions arrays when the charging energy becomes 
large. The nature of our phase corresponds to when lo- 
calization of single atoms in the wells destroys the super- 
fluid order parameter. Our transition is necessarily first 
order as the superfluid order parameter, and hence the 
superfluid response, jumps discontinuously at the transi- 
tion. The reason our transition necessarily destroys the 
superfluid order parameter is that, since our transition 
happens when the coherence length is larger than the 
well size or the interwell separation, boson localization 
introduces a length scale, the lattice constant, which is 
smaller than the coherence length for our transition. This 
leads to an energy scale larger than the superfluid order- 
ing scale. 

A Mott phase transition has recently been observed 
experimentally but the nature of the transition and 
of the Mott phase were not determined from such ex- 
periments. On the other hand, recent experiments arc 
consistent with a discontinuous transition irrespective of 
commensurability |^. 



In the present article we present the details leading to 
the conclusions of our earlier work. We present for the 
first time the explicit calculations, and we expand on the 
nature of the physics that controls when our transition 
is expected to happen as opposed to the transition pre- 
dicted in the early work |J| is expected to take place. 
Since issues of superfluid ordering and the different en- 
ergy scales in Bose Einstein Condensates (BECs) are es- 
sential for a microscopic understanding we start the arti- 
cle with sections that review well known physics of BECs 
and superfluids and their different length scales. We then 
study the nature of tunneling in optical lattices to review 
old results of ours 3] on the nature of tunneling in lat- 
tices of bosonic systems. When Josephson tunneling is 
suppressed before superfluid correlations are destroyed, 
the well known transition |2j follows. When superfluid 
correlations are destroyed before Josephson tunneling is 
suppressed, the transition we recently discovered follows 
. We then move on to estimate the properties of the 
Mott phase and the energetics of when the Mott phase 
wins over the superfluid phase leading to a localization 
transition. 



II. SIMPLE BOSONIC FLUID SYSTEMS 

In the present section we will review for completeness 
well known results in BECs and superfluid bosonic sys- 
tems. Most of this knowledge is due to studies of su- 
perfluid He'* which predates the experimental realiza- 
tion of artificially engineered BECs 7, 9} ■ These are 
superfluid systems as they posses a finite sound speed 
E0> ^3 suppressed long-wavelength scattering 
|13l The superfluid phase has off-diagonal order 

with breaking of gauge invariance [H, i.e. U{1) phase 
invariance, characterized by a coherent ground state of 
Bogolyubov pairs . 

Even though much of the material in this section is 
far from new, we included it in the article for complete- 
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ness and perspective for the subsequent sections. We also 
hope to give a flavor and emphasize the ideas of macro- 
scopic exactness that characterize stable quantum and 
thermodynamic phases of matter. While this perspective 
is not new, it is hard to tell if its importance is properly 
recognized in the BEC and atomic physics community. 
It is our wish to stress its importance. 



A. Noninteracting Boson Gas 

A gas of noninteracting bosons is described by the 
Hamiltonian 



(1) 



where the operators and at follow the harmonic oscil- 
lator commutator relations 

[at, at] = [aj:, ap] = ; [a^, at] = 5^^ (2) 

In the continuum limit, the sum becomes and integral 
and eo = 0. The dispersion relation is 



_ ^ 
2m 



(3) 



While in the continuum limit \kmin = 0|, in a finite box 
of size L with periodic boundary conditions 



k=—x+—y 



(4) 



where Hx, riy, and are integers different from 0. Thus 

in the box |A?mm| = ^^-^ and there is a gap that vanishes 
in the thermodynamic limit, that is, it vanishes as the 
number of bosons goes to infinity such that the density is 
fixed. For an optical lattice system the boson dispersion 
at long wavelengths is exactly as the one above for an 
appropriately chosen mass. 

The ground state wavefunction is trivial. 



|0) 



(5) 



where ao = . . This wavefunction refiects the fact 
that bosons can occupy the same state, so the ground 
state will consists of all N particles in the minimum en- 
ergy state. This state has energy eo = in the continuum 
limit. All the particles are coherent and have the "same 
phase": Bose Einstein condensation has occurred. This 
is the preamble to the occurrence of the phenomenon of 
superfiuidity. Bose condensation is a necessary condition 
for superfluidity to occur, but it is not a sufficient condi- 
tion. Bose condensation implies that, when the contin- 
uum limit is an appropriate approximation, the system 
is a quantum fluid. As explained below, a noninteract- 
ing Bose condensate lacks the necessary rigidity of the 
ground state wavefunction. A rigid ground state wave- 
function means that the dispersion relation is such that 



as k ^ 0, the excitation of quasiparticles costs an energy 
higher than Dissipation processes are prohibited 

when the system has enough rigidity as they would vio- 
late energy momentum conservation. Dissipationlcss re- 
sponse is the definition of superfiuidity as the coherence 
of the condensate can then survive weak enough pertur- 
bations. This necessary rigidity is provided by repulsive 
interactions among the bosons comprising the superfluid. 
In the absence of such interactions, the lack of rigidity 
permits scattering processes to occur which degrade co- 
herence no matter how weakly the system is perturbed. 



B. Interacting Boson Gcis 

BECs are really systems of bosons with repulsive in- 
teractions among them. They are supersaturated quan- 
tum vapors because they are metastable. They behave 
like quantum mechanical ground states for some time, 
but due to interactions or trap eff'ects they loose stability 
and decay. Since they live long enough to do measure- 
ments of their properties, and we can do calculations that 
agree with experiments as if they were ground states, it 
is not incorrect to treat them as ground states for times 
longer than their lifetimes. The finite lifetimes are both 
due to increasing interactions and trap effects. Usually 
their lifetimes are "long enough" that a lot of "ground 
state" physics can be studied. We concentrate in this 
latter physics. 

We will look here at some of the properties of BECs. 
Another system that can be described in an exactly the 
same manner and is actually a ground state is superfluid 
He^. The fact that is a liquid is not contradictory since 
the elementary excitations of the liquid have to repel each 
other as a stability condition. Otherwise, the elementary 
excitations can condense leading to some other phase of 
matter. 

We now proceed to describe the ground state and low 
energy excitation properties of the interacting Bose gas. 
We will concentrate on a boson system on a lattice with 
onsite Hubbard repulsion as we have bosons systems 
in optical lattices in mind. We do emphasize that the 
physics we describe here is universal as long as the sys- 
tem is a quantum fluid, i.e. there exists a macroscopically 
occupied lowest momentum state, and there is repulsion 
among the long wavelength elementary excitations of the 
system. When these two conditions are met, the physics 
is independent of the microscopic details of your system. 

We now concentrate on the Bose Hubbard Hamiltonian 
in a lattice (which of course can be an optical lattice) 

-t 



<13> 



+ U^ala^{a\a,-l)-e^a\ai (6) 



where is a chemical potential, U is the repulsion be- 
tween bosons, t is the tunneling amplitude between sites 



3 



or hopping term, and s represents the depth of the po- 
tential well. In a Mott insulating phase with localized 
bosons /i is irrelevant and can be set equal to 0. It is 
needed to describe the superfluid phase, as the approxi- 
mations made in the Bogolyubov Hamiltonian make the 
number of particles to not be conserved. In that case fi 
is chosen in such a way that the average number of par- 
ticles equals the total number of particles in the system. 
In the thermodynamic limit, the fluctuations around the 
average number of particles become negligible compared 
to the average. We emphasize that the nonconservation 
of particles in the Bogolyubov approximation is an arti- 
fact of the approximation. Physically the total number 
of bosons in the system is, of course, constant, but the 
number of bosons in the condensate cannot be constant 
due to spontaneous breaking of phase invariance|l5j. 

We consider a system with M bosons and N lattice 
sites. The number of bosons need not be commensu- 
rate with the lattice, i.e. M/N need not be an integer. 
The Hubbard interaction term in the Hamiltonian is after 
Fourier transforming 

= «r or at a.r 

X S{q + ki - fcg) S{k2 -h-q) (7) 
Performing some momentum summations we get to 



Tir = — y^ y^ at ar at ar 



(8) 



kiM 1 



The kinetic energy is a sum over pairs of nearest neigh- 
bors which describes tunneling from site to site: 



5A 



5,i 



where S is one of the vectors from a particle to its nearest 
neighbor. Fourier transforming diagonalizes the kinetic 
energy giving 



e atar 



k k 



5,k 



= H ^k'^k 



where er is the dispersion 



-ik-rx 



(10) 



(11) 



for a cubic lattice, with the energy measured from the 
state of zero momentum 

eg = — 2i(cos k^a + cos kya + cos k^a) + 6t (12) 



This shift of the zero of energy is irrelevant to the physics 
as it can be absorbed in a redefinition of the chemical 
potential. The dispersion goes like ta^k'^ = h^k^ /2m in 
the long wavelength limit. We see that efRcient tunneling 
corresponds to light particles and suppressed tunneling to 
very heavy ones. The full interacting boson Hamiltonian 
is then 

U 



'H — — at ^ar at jir 

kiM 9 

+ E ^k'^ln 'U^ + £ + U)Y. at a- 

k k 



(13) 



Now that we have ended with the Bose Hubbard 
Hamiltonian, we will consider ground state and low en- 
ergy excitation properties. At low temperatures a macro- 
scopic number of bosons occupy the zero momentum 
state as long as U is not too high. This is what it means 
to be Bose condensed. We take the number of bosons 
in the condensate to be Mq. A quantum state with a 
large quantum number behaves classically. An operator 
acting on that state corresponding to the large quan- 
tum number can be treated as a c— number, for its fluc- 
tuations are negligible and its quantum behavior is not 
seen. It is then possible to replace qq and aj by V Mq , 
the square root of the total number of bosons in the con- 
densate. For completeness we review all the tedious alge- 
braic manipulations leading to the reduced Bogolyubov 
Hamiltonian p. IT3| in appendix El 



M2 



MqU 
N 



where = 



fe/o k^O 

(14) 

f 4:Mf)U/N — fi and e and U factors have 
been subsumed in /i as the fixing of the average particle 
number will cause them to be swallowed by the chemical 
potential. 

For the ground state, which is a minimum of energy, 
the expectation value of Ti. must satisfy 



d{n) 

dMn 







If this condition is not satisfied, then by connecting a 
reservoir to our system we would have a flow of particles 
into and out of the system until equilibrium is reached. 
As long as the low energy eigenstates do not deplete the 
condensate, which they do not for not too high U, they 
satisfy the same extremum condition. Taking the partial 
derivative we find 



and 



Ll = U 

^ N 



^k-^k + -J^ 



U 



(15) 



(16) 
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The Hamiltonian in equation (|14|l can be diagonalized 
by making a Bogolyubov transformation 

h = + «fc«Lf bl = u^ai + v^a_j: (17) 



where we require 



(18) 



in order for the transformation to be canonical. The 
diagonal Hamiltonian has the form 

k 

where C is a constant. Substituting IjlTf) in the above 
expression we find 

k 

+ Ej:vl) + C (19) 
Comparing this expression with (|14|) we identify 

Solving these two equations together with u| — w| = 1 
yields 




^fe - ^£ - 



(21) 



This is the quasiparticle excitation spectrum. As fc ^ 

(22) 



Ek 



AMqU 2MoU 
erAl — 7^ =P\ 



N 



measured to have finite sound speeds^3| and are thus 
superfiuid. 

The ground state wavefunction is worked out appendix 
iBland found to be 

|*o) = n -e-^'"'^'''^"^''-^ {alrm (24) 



The factor (aj)*^'' represents the bosons in the conden- 



sate. In the absence of interactions, = 1, u^j = 0, so 
that in exciting bosons, one changes atomic boson num- 
ber by one. In the presence of interactions, states with 
momentum k and —k are mixed into the ground state, 
so when an excitation with momentum k is created, an- 
other with momentum — fc is destroyed. The excitations 
are coherent superpositions of atoms and "holes" . This 
is due to the nature of the coupling between k and —k 
in the Bogolyubov operators. With interactions both u-^ 
and increase. 

The ground state contains nonseparable, nontrivial 
correlations between the bosons that make up the system, 
i. e., it is entangled. The correlations in the ground state 
have entangled the bosons into a coherent state for the 
lowest energy state. The entanglement is so extreme that 
the bosons that make up the system cannot be excited 
at long wavenumbers. Their existence at low energies is 
impossible. Only sound can be excited, i.e. the excita- 
tions are Bogolyubov quasiparticles which do not resem- 
ble bosons whatsoever at low energies. In this limit, the 
principles of quantum hydrodynamics become exact be- 
cause the system does not dissipate and sound excitations 
become exact eigenstates of the system. The boson fiuid 
responds only collectively at long wavelengths. Aspects 
of the universality of superfiuid physics are presented in 
appendix |0 

We now proceed to calculate the self-consistency con- 
dition. Using the ground state wave function 124|) the 
self-consistency condition is 



^(vl/o|4a,,|*o) (25) 

fe#0 



Nm 



Let \a) = arj^'o), that is 



where p — ^2me^ = ^/2mt ak — hk. Since as fc ^ 
0, Ej^ — pCs, i.e. low energy excitations are sound, we 
identify the speed of sound to be 



2MoU 
Nm 



(23) 



This is a finite measurable quantity that indicates the 
presence of correlations among the constituents of the 
material. Such correlations stabilize persistent currents 
and give rigidity to the ground state by making it energet- 
ically expensive to excite the fundamental bosons inde- 
pendently. Instead, the excitations are collective modes 
and the system superfiows. The ground state has ac- 
quired rigidity. Artificially engineered BECs have been 



nl ^— O,^ Q,^ 1 — Q,^ 

— e -"'—e " -HO) (26) 

k'y^k 

The commutation relation [ar,a-l = 1 implies ar = 
Using this result 



(27) 



Substituting the Bogolyubov transformation a' - 
urb^ - - vrbr. 



^k"k 



vr. 



- ^A^)l*o> - ~Vkbl-,\^o) (28) 
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Thus (a I a) = w| and the self-consistency condition be- 



comes 



Er. 



k k 



nT.[- -l)+Mo (29) 



Ek 



This is the "particle conservation sum rule" . 

We will now make estimates of the number of particles 
in the condensate as a function of U/ 1. We use the result 

m 



1 



^k-2 



eg + 2{Mo/N)U 



/^,^ef: + 4{Mo/N)U 
Similarly we will also estimate 



- 1 



-I 



eg + 2{Mo/N)U 
/iI^eg + 4(Afo/iV);7 



(30) 



(31) 



We will concentrate on the long wavelength behavior, 
where equation H12|l for the dispersion gives eg ~ tk^a?. 

In this limit, the sum over k will be approximated with 
an integral: Y.k ^ {Na^/2^iT^) J dk dfl. Here N 
is the number of lattice sites. Notice further that for 
eg ~ tk^a^ < 2UMo/N, wg 7^ and ug ^ 1. On the 
other hand, for ifc^a^ >> 2UMo/N, wg ^ and wg ~ 1. 
This means that the system behaves like a superfluid for 
long wavelengths and like a free Bose gas for short wave- 
lengths. Superfluidity is thus a low energy phenomenon. 
This behavior can be approximated with the introduc- 
tion of a cutoff, so that above the cutoff wg ~ and 
Mg ~ 1 while below the cutoff they are given by the inte- 
gral through the region < k < kc- The cutoff is given 

by 



tk^a 



2^2 



,UMo 
' N 



V2 [UMq 

'W 



(32) 



For k < kc the system responds like a superfluid while 
otherwise it behaves like a dissipative free Bose gas. The 
important detail is that as long as the cutoff wavevector is 
smaller than the absolute cutoff from the lattice, ir/a, so 
that < 1, it will define the cutoff for all momentum 

integrals is fcc- The system generated a cutoff scale from 
the physics that lead to superfluidity. This scale is related 
to the all important coherence length of the superfluid 
about which we will have a lot to say later. Going back 
to ijSni and writing eg^ for 2UMq/N 




(33) 



and 




In terms of the cutoff 



-I 



Thus 



lQ^/2^T^ 



V2 k 



-I 



8^/2^l' 



-Na-'k 



31.3 



(34) 



(35) 



(36) 



Substituting this into (|29|l gives for the self-consistency 
condition 



M 



Na^k^ + Mo 



(37) 



8\/27r2 

This provides an expression for U/t in terms of Mq/M: 



1 = 



1 



8V27r2 
1 



U 

T 



8N/27r2 

(8V2^2)2/3 



" Af M 
2MoU\^''^ N Mq 



Nt 
N 



2Mo V N 



M 



2/3 



M 



1 



Mo 
M 



2/3 



(38) 



We now proceed to calculate the expected value of the 
Bose Hubbard Hamiltonian in the superfluid ground state 
in order to compare with the energy of the other phases 
and see which one is energetically more favorable. One 
of the main concerns of the present work is the Mott 
transition. So in particular we will compare this ground 
state energy of the superfluid phase with a certain Mott 
phase not previously studied in the literature^. The 
full Bose Hubbard Hamiltonian is given in equation l|13|l . 
The expectation value of the energy in the superfluid 
ground state is then 



(H) = -(e + ;7)M + ^egi;| 



U 



— V V(at a.- at ag ) - 6tM (39) 



ki ,fe2 



Notice that ar ar 

ki k2—q 



6r r - 

ki,k2—q 



at _ap . Using this 

k2-q ^ 



result we obtain for the interaction term 
i'Hi) = TT at ar + — y^(at Jit Jit Of 

^ -'' TV ^ k2 N ^ ^ ki+q k2-q '=1 ^2 



k2,q 



fel ,fe2 



U 
N 



E'S^ (at at ar ar 
/-^^ ki+q k2-q "^l '^a 



(40) 



fel ,fe2 
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On the other hand we have from (jDl|l Joining the two results we finally get 

%%l^o) = ^(1 + ^)'5s„_sJ*o) 

Similarly 

— v'i ut , - 

fci+q fc2-g ^ ^ Mr - I /ti, /i,2 

Kl+I? K2 — 9 

+ -Wr , -(^ol r r 1 - 



(at ^at ^ar ar ) 



fc2 fci+g 
K2 ki+q 



(1 



k-2 fci+g 



(1 



ki K2 



w2 It- 



)(1 



)(1 



k2-q 
v'i _ur I - 

k2-q ki+q 

ut - 
K2-9 



)5r r + vr vr vr -vr , -(^ 7"* i- , -^^ - 6^ - 

' fel,-fc2 fel fe2 fe2-? ki+q\ -ki-q -k2+q -ki -k2 



)'^fcl,-fc2 



^ki'^k2^k2-q'"ki+q^^ -k2+q-ki^ -ki-q,-k2 ^ '^-£2+'?, -fe2 "^-fel -g,-fel ) 



r 



(41) 



Thus performing the sums 

Ey^(at at ar ar 
^-^^ ki+q k2~q "^l ' 

^k 



= E-— 



k,ki 



(42) 



where ki = k — q. From H33|) and (|34|l = 1. Hence 



4^ Mj 23^3 

k 



d^k = —Na^kl 



(43) 



Usage of (1211) and (g^l into lHU gives 

E\^(at at ar ar ) 
/-^^ fci+g fc2-g fci ^2' 

fci,fc2 9 

~ U ^ 671^ 32/ 4^4 
As for the kinetic energy, with eg ~ tk^a^ and (|33|l 

> erwr — ^ — fcc / fc A: 

^ 57,5 



(44) 



16V27I 



(45) 



Substituting JSJ, gDJ) and gHJ into gives 

1 



■eM- 



16%/27r' 



-iVto^'fc; 



5^5 



NU - H ^ H — f 

\9 6V2 32; 4^4 



- 6tM 
1 



1 



16V27r2 
1 1 



\ Nt 

NU /2MoU 
9 ' 6V2 ' 32; 4^ Nt 
6tM 



(46) 



If the number of bosons is an integer multiple of the num- 
ber of lattice sites, the system is said to be commensu- 
rate. On the other hand, if the number of bosons is not 
an integer multiple of the number of lattice sites, the sys- 
tem is said to be incommensurate. Notice that the above 
result is general and thus applies for both incommensu- 
rate and commensurate systems. 

In the limit of small U/t, Mq ~ M and the second 
and third terms are small. The energetically expensive 
term, UN is small. There is an energy gain in becom- 
ing superfluid. As U/t increases ^^{^ goes to a constant 
as we will see in a succeeding section, and the second 
and third become of the order of A. If the second and 
third terms cancel the kinetic energy term —6tM, the 
system is not a quantum fluid anymore. It will thus have 
a transition into an insulating phase, a Mott phase. The 
commensurability does not play a role in the depletion of 
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the condensate, contrary to what was found in the earher 
work (21 . The earher work is not incorrect. In that work, 
the Mott transition happens when the superfluid within 
each weh dephases, preventing Josephson tunnehng and 
leading to a Mott-phase. That transition is continuous 
and only possible for commensurate values. For this to 
happen there needs to exist a superfluid within each well, 
which requires enough bosons per well (broken symme- 
tries can only happen for a macroscopic number of parti- 
cles) and the coherence length to be longer than the well. 
So in certain physical situations there is a Mott transition 
of a different type, discontinuous and incommensurate, 
since the superfluid does not survive within each well, 
and individual bosons are localized in each well. We will 
elucidate these physics in following sections. 



III. ORDER PARAMETER IN THE 
SUPERFLUID 

In the present section we go on with the study of the 
universal properties of the superfluid. We concentrate in 
what is the superfluid order parameter. This section also 
provides a connection with the new kind of Mott transi- 
tion we have uncovered in our work. In particular, since 
we will find that the superfluid order parameter increases 
as we increase the onsite repulsion U, i.e. as we approach 
a Mott transition, a transition in which the superfluid or- 
der parameter is destroyed must be discontinuous. The 
unique experimental signature of this is a discontinuous 
jump in any measure of the superfluid response. 

We have seen that the well depth e is subsumed into 
the chemical potential upon minimization of the ground 
state energy with respect to the number of particles in the 
condensate. The physical reason for this is that whether 
the lattice is deep or shallow is irrelevant as long as the 
system is superfluid; it superflows through it. Nonethe- 
less, the well depth is relevant to the Mott transitions and 
phases. Also, the kinetic term in ((BJ depends on e. This 
is so because decreasing the depth of the potential well, 
i.e. increasing e, will obviously decrease t as the barrier 
through which tunneling must occur is higher. The only 
vestige of e remaining in the superfluid is through t{e). 
We also notice that U/t oo with decreasing t. This 
however does not kill the order parameter of the BEC as 
long as the condensate is not depleted or Josephson tun- 
neling suppressed. If the condensate is depleted there will 
be a Mott transition in the universality class of the one 
uncovered in the present work. If Josephson tunneling is 
suppressed without destruction of the order parameter, 
there will be a Mott transition in the universality class 
discovered in the previous class ,2j . 

It is commonly said in the literature that the existence 
of a condensate is essential for superfluidity. The macro- 
scopic occupation of the lowest momentum state, i.e. the 
presence of the condensate, makes the boson system into 
a quantum fluid as long as the separation between energy 
levels is unmeasurable. The existence of the condensate 



does not imply the existence of a superfluid. Repulsive 
interactions imply the existence of superfluidity as long as 
there is a condensate, i.e. the system is fluid. We have 
seen that Bogolyubov correlations are the consequence 
of the interactions which leads to the necessary rigidity 
for superfluidity. Therefore, the correct order parameter 
that characterizes a superfluid is (a^a ^j), which is non 
zero only when there are Bogolyubov correlations. 

Now that we have recognized the correct order param- 
eter, the direct relationship with BCS superconductivity 
is obvious. The superfluid state is universal. Be it be- 
cause of a BEC that acquires rigidity due to the presence 
of repulsive interactions, or be it because of a Fermi gas 
acquiring rigidity due to attractive interactions, once the 
system orders into a superfluid, the macroscopic end re- 
sult in both systems is the same, independent of micro- 
scopies, modulo charging effects. There is an energy gain 
by ordering which stabilizes superflow and in both cases 
it is caused by the action of interactions. Both systems 
exhibit a Meissner effect: the BCS ground state does not 
respond to small enough magnetic flelds while the un- 
charged superfluid BEC system does not respond to slow 
enough rotations. 
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FIG. 1: (a) Dependence of the order parameter on the frac- 
tion of bosons in the condensate, (b) U/t vs. Mq/M 

The normalized Bogolyubov superfluid ground state 
wavefunction is 

|*o) = n -e-^'''^'''^^''-'^ {air" |0) (47) 



Thus the expectation value of the superfluid order pa- 
rameter is given by 



1 



where 



(48) 



v^(^)"(4a\-r 



(49) 
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Since [a_^aj:, (at)"(a^-)"] = n(aI.a^-)"-i we find 
1 



,t„t \n-l 



since Er 

k 



2 \"l 

k 



n. 



{n-l)\ 

(50) 



, 4_ 2Moa 



2MoU 
N 



(52) 



(51) 



or with ~ tPa"^ and 



Afo/M)2/3 



fc2a2 + 


/ 8\/27r2Af' 
V ^ . 












fc2a2 + 


8\/27r2Af' 




- Mo/M f/^ + ka\jk'^a? - 


F2(«^-^*^' 




Mo/Af)2/3 



As the potential depth is increased, the tunnehng am- 
phtude, and hence t in the Hamihonian, decreases. This 
does not degrade the order parameter. We also expect 
this to stabilize energetically the Mott insulating phase 
in order to minimize local number fluctuations. We thus 
expect a discontinuous transition into a Mott insulating 
phase, where the order parameter experiences a sudden 
jump to and the individual bosons, not the pairs, get 
localized into the potential wells. 

A superfluid has a density of Bogolyubov pairs. This 
means that even if only a few bosons are paired, the 
whole system superflows because of the rigidity acquired 
by pairing. We plot the dependence of the order param- 
eter on the number of bosons in the condensate and on 
U/t on figure n 



IV. MOTT INSULATOR AND MEAN FIELD 
SOLUTIONS 

The superfluid-Mott transition is between a superfluid 
and a quantum solid. As we have mentioned before in 
the present article and in earlier work 1] there are two 
types of Mott phases and Mott transitions. One of them 
exists when the superfluid survives within each well, both 
Josephson tunneling is suppressed within each well lead- 
ing to insulating behavior. The other which we have 
uncovered and it is partly the purpose of this work to 



explain occurs when the onsite repulsion destroys the su- 
perfluid order parameter and leads to boson localization 
in the wells. In systems without disorder, fluid-solid tran- 
sitions are usually discontinuous as the order parameters 
are too different 21\. In the superfluid system that is 
certainly the case. Since our Mott insulating phase is 
characterized by a well defined number of bosons per site, 
one could take the density as its order parameter. As it is 
very hard to have continuous transitions between a solid 
and a liquid, the continuous transition studied in the ear- 
lier work 13] only happens for a commensurate number 
of bosons per site. This is true, but whether a discontin- 
uous transition can also occur at incommensurate values 
is not so clear. In a following section we will study tun- 
neling between wells in order to elucidate the difference 
of the two Mott phases and transitions and speculate as 
to whether there might be an incommensurate transition 
in the case of suppression of Josephson tunneling. 

The Hamiltonian ^ with strong enough repulsion has 
a Mott insulating ground state. We rewrite equation ® 
as 

(5,i i i 

In the case where M < N there is no site with dou- 
ble occupancy. States with double occupancy are excited 
states. This is not true if M > TV. We treat the two 
cases separately. The strong repulsion ground state, or 
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Mott insulator, wavefunction corresponds to a rigid lat- 
tice with certain number of boson per lattice site, de- 
pending on the incommensuration. For M < it is 

M 

IV'o) = n4lO) (55) 

The double subscript notation helps differentiate the sites 
with bosons present (which are filled as j runs from 1 to 
M) from the empty sites (the ones left after all the values 
of j have been covered). The ground state wavefunction 
for pN < M < {p+ 1)N is 

AI-pN N 

l'/'o>= it 4.n(«^)'|0) (56) 
j=i 1=1 

With the Hamiltonian given by (I54|) the energy of the 
Mott state is found to be 

-{e + U) ^(V'ol^.lV'o) + ;7^(V^o|«.'|^o) (57) 

i i 

The first term corresponds to the kinetic energy. It de- 
stroys a boson in site i and creates it in its nearest neigh- 
bor site i + 6. Thus the overlap of the state resulting 
from this operation in the ground state with the original 
ground state gives 0. For M < N usage of the ground 
state wavefunction (|55|l in the above equation gives 

(V'olHlV'o) =-{£ + U)M + UM = sM (58) 

On the other hand, for pN < M < {p + l)N, usage of 
the ground state wavefunction (|56|) gives 

{^o\n\4>o) = - is + U)M + UiM - pN){p + i f 
+ U[N - [M - pN)]p^ 



(V'olHlV'o) = ~£M + 2pUM -{p+ l)pNU (59) 

This energy is to be compared with the ground state en- 
ergy of the superfluid H46II . For U/t small, we see that 
the interaction terms in the superfluid are almost irrele- 
vant when compared with those of the Mott phase so the 
superfluid state is the stable one. We also notice that for 
not too large an M the interaction terms in the superfluid 
grow a lot faster than the interaction terms in the Mott 
insulating phase. The physical reason for this behavior is 
that in the fluid all particles overlap with each other as 
there is a probability for each of them to be at each lattice 
site. Therefore for some critical U/t the Mott insulating 
phase will become stable causing the system to undergo 
a quantum phase transition. Whether M is commen- 
surate or incommensurate with the lattice is irrelevant 



to the energetics: The state with all bosons localized is 
energetically favorable to the state with all bosons super- 
fluid. Whether the instability of the superfluid is toward 
a phase with all bosons localized, or into a state with 
coexistence of a commensurate Mott insulator with the 
leftover bosons superfluid is not straightforward to an- 
swer variationally. It can be answered experimentally 
as in the case of coexistence the superfluid response will 
have a discontinuous jump to a smaller nonzero value, 
while in the case of all bosons localized the superfluid 
response will have a discontinuous jump to 0. 

Let us take a look at the expectation value of the order 
parameter in the Mott phase. It is given by 

(«£«-£) = ^ E e*('-'-^^-)(V'o|a.a,|V'o) (60) 

<i]> 

where 

{Q\X{afa,a,\{{a]r\Q)^ \{{m,af{a\ym (61) 

In the expression above, qi indicates the number of 
bosons in site I and Fij is a function of bosonic oper- 
ators on sites i and j . 

n [mAiAm = n {mA''m{a\r\o) 

+ \{{Q\F,,ar\a\ram 

(62) 

where we have used ai — . Continuing the process we 
find 

= Hqmafafa^aMn^y'lO) (63) 

Since there is a non matching number of operators acting 
on the ground state this yields 0. The order parameter 
is in the Mott phase 

= (64) 

Actually this is obvious because the Mott state has a well 
defined number of bosons per site and number phase un- 
certainty makes the off-diagonal order zero. Since as U 
increases, the order parameter increases in the superfluid 
phase and the order parameter is zero in the Mott phase, 
the only way to go into the Mott insulating phase is to 
have a discontinuous change in the superfluid order pa- 
rameter. The Mott transition is discontinuous, i.e. first 
order, with the superfluid order parameter, and thus the 
superfluid response, experiencing a sudden jump to 0. 
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V. INCOMMENSURATE MOTT PHASE 

The problem of the phase transition in a BEC has been 
treated in the not distant past by M.P. Fisher, et. al. 0- 
Their treatment rehes in having a Bose condensate within 
each well and the order parameter (a^a-k) being fixed 
while the phase fluctuates more and more as the super- 
fluid density stiffens up due to the increase in U /t. It is 
interesting to note that despite the increase of the order 
parameter with decreasing t, the speed of sound will be 
collapsing. This is the loss of rigidity due to phase fluc- 
tuations that we just mentioned, which signals a transi- 
tion into a Mott phase different than the Mott insulat- 
ing phase studied in the last section. The Mott phase of 
Fisher, et. al. will be a solid with no long range coherence 
but superfluid coherence within each well. The transition 
out of the superfluid could be a continuous transition as 
the phases at different lattice sites lose coherence and cor- 
responds to the transition studied by Fisher et. al.. It 
might be possible for this transition to be discontinuous 
for incommensurate number of bosons per lattice sites. 
We will try to study this matter in the near future. 

While this transition can certainly occur, the transi- 
tion which takes place in most experiments in optical 
lattices corresponds to individual boson localization into 
a Mott insulating phase as studied in the previous sec- 
tion. The reason for this is that most experiments have 
low superfluid densities and equivalently, small number 
of bosons per site 0, 0] . Since one needs a macroscopic 
number of bosons to have a superfluid and one needs 
coherence lengths shorter than the well size in order for 
superfluidity to survive within each well. The experimen- 
tal conditions of most lattice experiments cannot support 
a transition as studied in the work of Fisher et. al. in 
which the superfluid order parameter does not die. 

In our transition, when the individual bosons localize, 
the Bogolyubov boson pair correlations are broken apart 
within the lattice and within each well. In a commensu- 
rate system, if the depth of the potential well is enough 
to overcome the on site repulsion J7, the system will lo- 
calize and each site will be filled with the same number of 
bosons. The maximum lattice vector will be — . For an 

a 

incommensurate system the situation is somewhat differ- 
ent. Suppose there are N lattice sites, but iV-|- 1 bosons. 
Then as the depth of the potential well is increased U 
increases and N bosons will be localized one per site, 
forming an underlying lattice. There are several possibil- 
ities for the remaining | bosons to go. Two of them are 
shown in the Figure |21 

If only onsite interactions are taken into account, the 
two configurations are degenerate. In the case of m extra 



bosons, the degeneracy is 



The interbosonic 



m\{N-m)\ ■ 

distance is different for each degenerate configuration. 
Since there is no preferred distance, the extra bosons 
will not be initially trapped. If the well depth keeps 
increasing, the extra bosons will eventually be trapped 
too. But they will do so with many different wavevectors 
(fi. The structure formed will be an amorphous solid. 



• V • 

• • • • 

• • • • 

Configuration T 




Configuration II 



FIG. 2: Two possible configurations for extra bosons in an 
incommensurate lattice. 



In real life, besides on site interactions, there are inter- 
actions acting at long distances. They can be attractive 
or repulsive. These interactions are weaker than the on 
site interaction J7, hence not important for a commen- 
surate system. But if there is incommensuration, these 
interactions will probably play a role in determining the 
positions of the extra bosons. Another possibility is that 
when increasing U the extra bosons get trapped in one 
of the possible configurations. If the configurations are 
metastable, then the system will remain there for a long 
time. 

In our previous work, we had determined that a Mott 
localization transition could happen for incommensurate 
values of the number of bosons per site . The nature 
of what this phase was was not elucidated. The diffi- 
culty lied in where to place the extra bosons after the 
integer part was subtracted from the total. As far as the 
Mott phase was concerned, they could fall in the wells 
at random as long as not two of them lie in the same 
well, for this would cost extra C/'s of energy. In real life, 
an extra interaction might cause them to order in a cer- 
tain way and decide the positions where they lie. We 
do not consider such a case here. Another possibility 
in the Bose-Hubbard model is that the integer part of 
the number of bosons per site (which is commensurate 
by definition) gets localized in each well and the leftover 
fraction constitute a Bose fluid. This would correspond 
to Mott-superfluid coexistence. We find that no such co- 
existence is possible. 

When the total number of bosons is very near a com- 
mensurate number, the few leftover bosons are not lo- 
calized and gain kinetic energy by being spread through 
the lattice. On the other hand, we do not think this 
constitutes a phase where we have coexistence of two 
phases. This is because the fraction of bosons which can 
be spread over the lattice collapses to zero as the lattice 
is made larger and larger. A true thermodynamic phase 
survives the thermodynamic limit. In fact, its properties 
become exact in the thermodynamic limit. Hence, the re- 
gion where the fraction of bosons is free over the lattice 
exists only due to finite size effects. We thus conclude 
that for incommensurate values there is a Mott phase 
where the fraction of bosons that causes the incommen- 
suration gets localized at random places in the lattice, 
forming a kind of amorphous solid. 

We now estimate how this comes about. We will study. 



11 



making mean field estimations, whether it is energetically 
favorable for the ground state to have all the bosons local- 
ized, or to have a commensurate number localized with 
the leftover fraction forming a gas or fluid coexisting with 
the Mott phase. The coexisting fluid will of course be a 
superfluid because of the repulsion, but we will only con- 
sider the fluid to be a Bose condensed free boson gas as a 
precondition for superfluidity is to have a Bose conden- 
sate. 

In both the Mott phase with all the bosons localized 
and the phase with coexistence of the BEC and the Mott 
insulator (MI), we have J^ii^i) = Therefore, the 
chemical potential term will not favor any of the phases. 
In this section it will be dropped from our Bose Hubbard 
Hamiltonian for energetic comparisons: 



<ij> 



(65) 



Without loss of generality we will consider the case where 
M lies between N and 2 A'', which we parametrize as M = 
{l+p)N with < p < 1. The wavefunction that describes 
the Mott phase with all particles localized is given by 



M 



i*M)=n4io) 



(66) 



where the jVs run over lattice sites and the lattice sites 
are filled such that there is not less than one and not 
more than two particles per site. This way we have the 
minimum Hubbard energy possible. The expected value 
of the kinetic energy in this state is zero. Hence, the 
expected energy of this state is 



GS 



(3M - 2N)U = UN{1 + 3p) 



(67) 



This follows straightforwardly because we have M — N 

sites with two particles at an energy cost AU , and 2N—M 
sites with one particle at an energy cost U. 

For the phase consisting of the coexistence of a com- 
mensurate Mott phase (iV particles with one particle per 
well) with a BEC of the leftover M — N particles the 
wavefunction is 



N 



|*MB) = (aS)^-^n«llO) 



(68) 



where is the operator that creates a boson with mo- 
mentum k 



0. al 



i'^/VN)Y.^al■ Only the fluid 
particles contribute to the expected value of the kinetic 
energy, which is given by 



{Eh) = -6t{M - N) 



(69) 



In order to estimate the expected value of the interaction 

energy, we need to make a few estimates. The Mott com- 
ponent of this coexistence phase has exactly one particle 



per well. The fluid component, of course, has a fluc- 
tuating number of particles per well, whose average is 
(M-iV)/A/' leading to 



(rii) = 1 + 



M-N M 



N 



N 



(70) 



In order to estimate the expected value of nf , we reckon 
the following way: the Mott component of the coexis- 
tence phase does not have particle number fluctuations 
per site, leading to (An?)M = 0. The BEC component 
of the coexistence phase has fluctuations per site lead- 
ing to typical fluctuations of a random liquid (classical 
or quantum) (An?)MB = (M - Nf/N. This number 
comes aboiit because there is a quantum amplitude for 
all the fluid particles to be at each well (this is what it 
means to be fluid or delocaHzed). Hence, all the parti- 
cles contribute to Anf such that Anf oc (M — A^)^. This 
needs to be multiplied by the probability of flnding a par- 
ticle in the well, which is The total fluctuations are 
given by quadrature: Anf = {Anf)M + {Anf)MB- From 
this we immediately obtain 



(71) 



N m 

The energy of the coexistence phase is 

E^i = -6tNp + UN[{N + 1)/ + 2p+l]. (72) 



We define the quantity 



pMB pM 
^GS ~ ^GS 

NU 



6t 
U 



+ l]p+{N + l)p' 



(73) 



which when positive, the Mott phase is stable and when 
negative, the coexistence phase is stable. It follows that 
for p between and {{6t/U) + 1)/{N+1) the coexistence 
phase is stable while for p > {{6t/U) + 1)/{N + 1) the 
incommensurate Mott phase is stable. 

We see that when N bosons get localized into a Mott 
phase with a single boson per site, the extra bosons one 
adds on top of these will gain more kinetic energy by de- 
localizing than the potential energy penalty coming from 
the overlap with the localized particles. Hence, they form 
a sort of fluid. The potential energy cost becomes too 
much when there is a fraction p = {{6t/U) + 1)/{N + 1). 
At this point they will localize too. This will happen 
despite there being an incommensurate number of them. 
In the absence of other long range interactions to tell 
the extra bosons where to localize, they will localize in 
random wells as long as no two of them go to the same 
well. For macroscopic lattices we see that the fraction 
of bosons that will delocalize becomes arbitrarily small. 
Therefore, this is a flnite lattice size effect and does not 
reflect the properties of a true thermodynamic phase of 
matter, which by deflnition becomes better defined in the 
macroscopic limit. 
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VI. JOSEPHSON VS. INDIVIDUAL BOSONS 
TUNNELING 

In order to clarify the difference between the Mott 
phase we have uncovered and that of earlier work, we 
study BEC tunneling in an optical lattice. It has been 
verified experimentally in one dimensional optical lattices 
that there is coherent quantum tunneling in BEC'sp^. 
We will try to take an in-depth look at tunneling in opti- 
cal lattices in a general number of dimensions. Of course, 
the more relevant cases to superfluid-Mott transitions are 
those of higher dimensions, as superfluid ordering cannot 
occur in one dimension (see Appendix . In order for 
the tunneling to be relevant, the hopping or tunneling en- 
ergy scale, t, needs to be larger than the kinetic energy 
scale for the whole lattice: 



Tl = 



(74) 



where m is the mass of the bosons and Q is the volume 
of the whole lattice. If the kinetic energy of the lattice 
overwhelms the hopping, the tight-binding approxima- 
tion with bosons hopping from lattice site to lattice site 
with tunneling t is no longer valid for delocalized bosons. 

On dimensional grounds one expects the tunneling rate 
to be proportional to t/h. On the other hand, the pro- 
portionality factor should be a universal dimensionless 
function of the number of lattice sites N , the number of 
bosons M and/or the number of bosons in the conden- 
sate Mq, and the different energy scales in the problem. 
Besides t and fi^ /2mn'^^^ , there is the on-site repulsion 
at each well U, the well depth e, and the kinetic energy 
scale of each well: 



T — 



2mF2/3 



(75) 



where V is the volume of each well ~ a'^ with a the 
interlattice spacing. We, of course, have NV = il. The 
wells are taken as large boxes. This approximation is 
valid as wells in optical lattices are fairly large (a ~ lOO's 
of nm) in typical lattices 

The interwell tunneling rate or particle current is 



t / t t T^, 



(76) 



The universal function F will behave differently in dif- 
ferent physical regimes. For example, if the wells are 
very deep, e ^ i, the tunneling rate should be expo- 
nentially suppressed t ~ exp(— ae/i) with a a positive 
numerical constant. We will never write the dependence 
of the function F on t/e because the only effect of e for 
the suppression of tunneling is through the decrease of 
t. The number of bosons in the condensate, Mq, is de- 
termined uniquely by the number of bosons M , and the 
repulsion U. We see that the relevant parameter that sets 
the scale for the interaction for the BEC is the combina- 
tion UMq/N. Hence the universal dimensionless function 



that determines the tunneling current enhancement does 
not depend on U, Mq and N independently, but rather 



f[Mq,nI,^^f(J^,^ 

' ^' ' U' U J \UMq' UMq 



(77) 



We remind the reader that the dispersion relation for 
the superfluid is 



4C/2M2 

7V2 




(78) 



where we make the continuum approximation e^j = 
h^k^ /2m. As /c — > the dispersion becomes phononic, 
which leads to dissipationless flow. The dispersion at 
large k becomes h^k"^ /2m, which leads to dissipative flow. 
Notice that the crossover betwee n dissipative a nd nondis- 
sipative flow happens at fc ~ ^/StvIJM^JWN . This de- 
fines the coherence length of the superfluid 



SmUM) 



(79) 



When the system is probed at length scales L > 
it responds collectively as a superfluid, while at length 
scales L < it behaves like a gas of independent bosonic 
particles, thus responding dissipatively, or even not be- 
having like a fluid at all if the spacing between levels is 
relevant. We note that an infinite system would become 
superfluid for any U no matter how small, as one can 
always go to a long enough length scale to see the system 
respond collectively as a superfluid. In here we do not 
have in mind infinite systems but finite wells of super- 
fluid. In this case of finite systems, if the kinetic energy 
scale of the well, T^, is greater than the effective inter- 
action, UMq/N, the system will not be a superfluid as 
the linear size of the system, V^^^, is less than the co- 
herence length of the superfluid. Independent of whether 
the bosons are superfluid in each well, the system could 
still be a superfluid within the whole lattice. 

For the full lattice, when the kinetic energy scale is set 
by tunneling between sites (e^: ~ tk'^a?), the coherence 
length is 



Nt 
AMqU 



(80) 



The size of the lattice is, of course, L — Na and the con- 
dition for superfluidity (^ < L) is U/t > I/AMqN. We 
see that either for a very large system (A^ ^ oo) or very 
large number of bosons (Mq — > oo), it becomes arbitrarily 
easy to become superfluid in the sense that an arbitrar- 
ily small U will order the boson fluid into a superfluid 
at low enough temperature. Conversely, as we make the 
system small, or we reduce the number of bosons, it be- 
comes harder to order into a superfluid, which will show 
up experimentally as increased dissipation rates. 

Even though a high enough U can order the system 
into a superfluid, too high a U will make it into a Mott 
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insulator 0, The transition to the Mott insulator 
can be continuous^ or discontinuous The continuous 
transition occurs when there is a superfluid0 in each well 
and U prevents Josephson tunneling between wells. The 
continuous transition only happens for a commensurate 
number of bosons per site. When the bosons in the wells 
are not superfluid within wells the transition is discon- 
tinuous and happens irrespective of whether the number 
of bosons per site is integer or notQ. The second transi- 
tion always happens when the number of bosons per site 
or U/t are sufficiently small as in these cases ^ > a. In 
fact from our estimate H8U|) there cannot be a superfluid 
within each well when 



Nt 
AMoU 



> 1 



(81) 



For this case the Mott transition will be like ours if this 
condition still holds at the critical value of U/t, other- 
wise the Mott transition will be that of suppression of 
Josephson tunneling 2] . 

We now proceed to consider the case when there is 
a superfluid within each well. This will elucidate the 
microscopic physics of how onsite repulsion suppresses 
Josephson tunneling, leading to the transition studied 
by Fisher, et. al. We will sec that commcnsuration 
plays no role in the suppression of tunneling. This leads 
us to suspect that a first order transition can occur for 
this case too contrary to the early work The reason 
this was not found out in that work is that they perform 
renormalization group studies which can only access con- 
tinuous transitions. Whether a discontinuous transition 
can occur in this case is a matter for further study. 

We evaluate the current operator (25i| 



J 



{N -N) 



dT 

2it y-^ 

k 



[N - N, Ht] 



(82) 



the number of atoms within the wells is large. For such a 
calculation the only relevant excited states are two quasi- 
particle states, one in each adjacent well, with opposite 



momenta, i.e. |n) = btb'' jtlV'o), with energies E„ 



This calculation was first done by B. D. Josephson for 
Cooper paired fermionic superconductors |25j . We will 
follow his calculation closely. When the wells are su- 
perfluid, the tunneling current is dominated by coherent 
tunneling of Bogolyubov correlated pairs. 

In order to calculate in linear response theory"^, we 
turn on Ht adiabatically from time T — —oo, starting 
from the ground state of the noninteracting system at 
such a time. If we expand the wavefunction in terms of 
the noninteracting eigenstates 

l*> = ""(^) HEnT/h)\n) 

n 

we obtain the well known formula 



a„(T) = 6n,o - ie:xjp{iEnT/h){n\Ht\ipo)/{iE„ 



(84) 



where the vacuum state energy Eq has been chosen to 
be 0. The parameter e was introduced to control the 
adiabatic switching of the interaction. The causal limit 
e ^ 0+ should be always understood. Our interaction is 
turned on fully at time T — and the wavefunction of 
the interacting system is given by |^'(0)) = \^o) + \6(p) 
where 



\6cp) 



E 



\n){n\Ht\%) 
En - ie 



(85) 



The Josephson current is given by (5'(0)| J|^(0)) = 
(0| J|(5(^)+ complex conjugate. With 



\s^) = tj2 



^bibK + bKbi 



(86) 



where 



Ht 



(83) 



is the tunneling part of the Hamiltonian and T denotes 
time. Now a^; and aj: are boson operators for different 
adjacent wells. At first we consider the tunneling current 
when we have a superfiuid within each well. We note that 
the ground state of the reduced Hamiltonian H14(l when 
t = is an outer product of the Bogolyubov ground states 
for each well, which we denote as our "vacuum" state 
\ipo). The excited states of such a system are denoted 
by \n) and are constructed by applying any number of 
Bogolyubov creation operators to the vacuum state. The 
energies of such states. En, will be the sum of the en- 
ergies of each of the operators. When t ^ there will 
be tunneling between adjacent wells, which can be calcu- 
lated accurately from linear response theory as long as 



where 9 is the phase difference between the superfluid 
order parameters in adjacent wells, we obtain 



J 



2*2 



•E 



9 2 
k k 



(87) 



analogous to the fermionic superconducting case [23 • We 
note that if 6* is there is no Josephson tunneling as the 
superfluids are in phase. Since u^Vj: = UMq/ Ej^N , we 
obtain for the enhancement factor F = Jh/t 



2tU^M§ 



smf 



E£;3 



In order to estimate this quantity we approximate 
the quasiparticle dispersion as phononic, E'i ~ 

2U MofP'k'^ / mN , and we cut off the high momenta at kc 
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of the order of the inverse coherence length. Converting 
the sum into an integral 



1 



1 fmNV^/^ 



^ El (27r)3 V 2C/Mo 



3/2 



(Pk 



(89) 



Note that this integral is infrared divergent but the di- 
vergence gets cut-off by the inverse linear well size V'^/^. 
We finally obtain for the enhancement 



F 



1 tN (UM,\^'\^AUM^ ^^^^ 



167r2 UMo V Ty,N 



T,„N 



In order for the system to be superfluid at all within a 
well we must have the linear dimensions of the well larger 
than the coherence length, y^/^ > ^, and there needs 
to be enough bosons in each well as to have superfluid 
broken symmetry. Otherwise, there will not be Joseph- 
son tunneling per se and the calculation of tunneling will 
have to be done for noninteracting bosons trapped in 
wells. This condition is UMq/T^N > 1/4 with smaller 
values corresponding to a physically impossible regime. 
When the conditions are such that there can be Joseph- 
son tunneling between sites, intrawell Coulomb repulsion 
will suppress Josephson tunneling, leading to a Mott in- 
sulating phase as studied originally ,2|. On the other 
hand, when Josephson tunneling cannot occur from well 
to well because the coherence length is longer than the 
well size, the Mott transition will be single boson local- 
ization as uncovered by us 0. The conditions in Mott 
transitions in optical lattice experiments 0, are such 
that the coherence length is too long and the number 
of bosons too few as to support Josephson tunneling be- 
tween adjacent wells as there would not be a superfluid 
within each well. 

We concentrate on the case when the linear size of the 
wells is too small {V^^^ < C)- this case the interaction 
cannot order the system into a superfluid within each 
well (although it can certainly be a superfluid within the 
whole lattice), and its only effect is a renormalization of 
the mass of the bosons. Hence the system is described by 
the Hamiltonian H83|l with m the renormalized mass of 
the bosons. There is no Josephson tunneling per se as the 
system is not superfluid. On the other hand, there can be 
coherent quantum tunneling of the (nonsupcrfluid) boson 
liquid between wells. 

The Hamiltonian (|83|l is diagonalized by the canonical 
operators 



^ _ H + H ^ _ H - H 



(91) 



to yield 

^ = E { (^S - + + i)A\h} (92) 

k 

The Heisenberg equations of motion lead to the following 
time dependent operators Aj^{T) = A^eyiY>[—i{ej:—t)T/li\ 



and A^{T) = A-j:eyi-p[—i{e^ + t)T /hi]. The common phase 
factor —ie^T/h will be dropped as only phase differences 
affect the physics. From these time dependent operators 
one immediately obtains the operators into which 
and their conjugates develop under Hamiltonian evolu- 
tion: 



a^(r) = COS 
a^(T) = cos 
aUT) = cos 



aJ(T) = cos 



tT 
tT 
tT 

tT\ _t . . 
— a- — z sm 
h J k 



a -~ — I sm 

k 



at 

k 



at 

k 



(93) 



The "barred" and "unbarred" operators are orthogonal 
raising and lowering operators as unitary Hamiltonian 
evolution is a canonical transformation thus preserving 
the commutation relations of the operators at time zero. 

We can now calculate the coherent tunneling current 
between two wells. We take the right well to have {N + 
M)/2 bosons and the left well with {N - M)/2 bosons. 
The initial wavcfunction is then 



l*(0)) = 



(„t)(iV+M)/2(gt)(iV-Af)/2 



{N/2 + M/2)\{N/2 ~ M/2)\ 



l^o) (94) 



where the subscript means the smallest momentum 
state into which the bosons condensed. The wavcfunc- 
tion at time T, |*(T)), has the exact same dependence 
on the operators a\{T), a\{T) that |5'(0)) has on aj, 
aj. In order to obtain the coherent tunneling current J 
we need to evaluate (AA^(T)) = (*(r)|AA^|*(T)) with 
AA^ = ajao — ajao counting the excess number of bosons 
in the right well. Using 

AiV = cos f — j (at(r)ao(T) - 4(T)ao(T)) 

+2zsin f — j (4(r)ao(T) - aS(T)ao(T)) (95) 



we obtain the tunneling current 
d 



J 



dT 



{AN{T)) 



2t ^ . f2tT\ 



(96) 



where the periodicity is an artifact of having a two site 
system. The proportionality of the tunneling current on 
the initial number difference of bosons between the sites, 
M, is not an artifact. This coherent quantum tunneling 
can and will be suppressed by strong enough repulsion 
within the wells, leading to boson localization within the 
wells. 

Before concluding we will make some order of mag- 
nitude estimate of when the Mott transition happens. 
Let us first consider the case when the coherence length 
is smaller than the well size and there is a macroscopic 
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number of bosons per well such that there can be a super- 
fluid within each well. In this case, the Mott transition 
will be suppression of Josephson tunneling as studied the- 
oretically before the advent of BECs "^l. In this case, the 
energy scale associated with tunneling is the Josephson 
tunneling energy 



1 t fUMo 



1/2 



In 



4UM, 







T,„N 



(97) 



where AUMq T^N > 1 for Josephson tunneling to be 
possible. When the average Hubbard repulsion energy 
per site, Ejj = U{Mo/N)'^ , is larger than the Josephson 
energy the Hubbard interaction prevents tunneling, mak- 
ing the system into a Mott insulator. Since we expect the 
logarithm not to deviate much from factors of order 1, 
we can estimate this condition to be about 



-1/2 



Mo 
N 



(98) 



For the case of long coherence length or few bosons per 
site, the transition will be individual boson localization 
with suppression of individual boson tunneling rather 
than Josephson tunneling. The energy scale for indi- 
vidual boson tunneling is set roughly by t. We would 
expect the transition to the Mott phase to happen when 
this energy is swamp by the Hubbard interaction energy 



U 
J 



N 
Mo 



(99) 



which is very roughly a condensate depletion condition. 
When the lattice tunneling is more relevant than the 
intra-well kinetic energy, t > Tw and the number of 
bosons per site is of order 1 , the boson localization phase 
transition described in this work will happen. On the 
other hand, iit< Tw or the number of bosons per site is 
considerably greater than one, one expects a Mott tran- 
sition with suppression of Josephson tunneling and su- 
perfluidity within each well. If the number of bosons per 
site is considerably less than one, one can probably not 
speak of either transition. 



is becoming small as U/t increases. As the condensate 
gets depleted the system is becoming "less fluid" and will 
solidify into a Mott insulating phase with a well defined 
number of bosons per site. When the coherence length 
is longer than the interwell spacing, or the number of 
bosons per well too small, the Bogolyubov order param- 
eter will be zero within each well and in the quantum 
solid phase. There will be a discontinuous jump in the 
superfluid response at the Mott transition. This transi- 
tion can happen regardless of commensuration and the 
conditions in most lattice experiments are conducive to 
this transition 

The physics of the transition we have described in the 
present work differs considerably from the one proposed 
in the pioneering theoretical work on Bose systems 
and thus requires some comments. The first and, per- 
haps, most important difference is that we do not con- 
sider a phase only model as it is not appropriate to the 
Mott insulator at low densities since a superfluid order 
parameter cannot survive within each well. In that orig- 
inal work, the phase of the superfluid was disordered by 
the increasing repulsion thus leading to a transition. In 
such a transition the solid would have a superfluid order 
parameter within each well, but the phase of the order 
parameter will have become randomized exactly analo- 
gous to what happens in Josephson junction arrays when 
the charging energy is sufficiently large. Such an insulat- 
ing phase does not correspond to the Mott insulator we 
studied here and cannot exist at small number of bosons 
per site as one needs a macroscopic number of particles 
to have a superfluid order parameter. It can also not ex- 
ist when the coherence length of the superfluid is longer 
than the well size. Finally, we stress that the phase only 
model studied in the early work 0, HO] is a correct model 
of an array of Josephson coupled superfluid systems and 
should work for bosons in an optical lattice when the 
number of bosons per site is large enough, and the co- 
herence length short enough as to make superfluid order 
within each well possible. Josephson coupled systems can 
easily be studied in an optical lattice [22, 11^ . It will be 
extremely interesting to see what the experimental phase 
diagram turns out to be. 



VII. CONCLUSION 



APPENDIX A: DERIVATION OF THE 
BOGOLYUBOV HAMILTONIAN 



In BECs in an optical lattice the Mott transition is in- 
duced by making the wells on the lattice sites deeper (e 
is the lattice depth) thus suppressing the tunelling am- 
plitude t between near neighbor sites, i.e. making the 
bosons "heavy" . Thus the relevant ratio of on-site repul- 
sion U to hopping t increases. As this ratio increases, 
the superfluid order parameter increases as shown in the 
figure^ The increase of the order parameter means that 
the superfluid phase is becoming more stable over the 
Bose condensed gas. On the other hand, we also plot 
the number of particles in the condensate. This number 



The interaction part of the Hamiltonian can be written 



k,k2 ,q 



k,k2 ,q 



E«b (Al) 

k,q 



16 



Now we start analyzing terms from 7i/. Separating from 
the momentum sums all the terms with oq or aj and re- 



placing the values of these operators by \/Mq since the 
condensate is macroscopically occupied and these opera- 
tors behave like c- numbers, we get 



I 



We concentrate on the universal ground state and low 
energy properties of the superfluid. We thus neglect the 
terms with factors of ^/Mq or smaller since the macro- 
scopic occupation of the condensate will make them 
unimportant. We arrive then at the reduced Bogolyubov 
Hamiltonian 

fc#0 

+ ^E(«R.+«-.^'^.^) (A3) 

where 

= e^: + (AUMo/N) - ^ . (A4) 

The extra term j;'^^ ^ataj: in equation I)A1|I can be 
written as U^j:ataj:. In this way, it is seen that such 
a term just renormalizes the chemical potential. It has 
therefore been absorbed into fj,. 

APPENDIX B: BOSONIC SUPERFLUID 
GROUND STATE 

The manipulations reviewed in the present appendix 
are original to Bogolyubov and we follow closely K. 
Huang 27] The ground state wavefunction for the Bo- 
golyubov Hamiltonian l|14|l has the form 

|*o) = X! E •■• (cni,miC„2,„i2 ... )|nimi; n2m2; ... ) 

rii ,mi n2 ,m2 

(Bl) 

where rii is the number of particles with momentum ki 
and rrii is the number of particles with momentum —ki. 
It is obviously the state with no excitations 

^a;I*o) = {u^aj: + z;,V_-)|*o) = 



Let Cn.m — Cni,miCn2,m2 ■■■ : ^ — '^^j ~ ^nd \n,m) — 
|nimi; n2m2', ... ). Then 

y^(M Cn,mVn \n-l,m) +V Cn^mVm + 1 1^1, TO -|- 1)) = 

(B2) 

Changing dummy variables n n + 1 in the first term 
of the sum and to — > m — 1 in the second term 

y^(M Cn+i,mVn + 1 + V Cn,m~iVrn) \n, to) ^ (B3) 

n,m 

For n ^ m and m = we have 

U Cn+l^rnVn + 1 + V C„^„i_i\/to = 

Cn+1,0 = for all ?i's 

for m = 1 

u Cn+i,iVn + 1 + V Cnfi = u c„+i4 sqrtn + 1 = 
c„+i,i = 

and so on for all Cn,m- Thus we only need to consider 
71 = TO in Cn,m- In that case we have from (jB3|l 

C7n,m ~t- V Cjn— l,m— 1 — 

_ V _ ( v 2 — r "^^3 

U U U 



c™,™ = (--)'" co,o (B4) 
u 

Substituting this into (|B1|) we find 



oo 
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where we can write the state composed of n particles with 
momentum k and n particles with momentum — fc as 



i4r i^Ur 



hence 



(B6) 



|0) (B7) 



|co,o| 



y (— ) 



k_\2n 



(B9) 



Remembering that X]n(l^')" ~ i~J^/u^ ~ ^\ arrive 
at the expression |co,oP — The normalized ground 
state wavefunction is then 



i*o)=n- 



(BIO) 



i*o)=n 



co,oe 



— ^^^^ ^ 



|0) 



(B8) 



which is a coherent state. We need to find the normal- 
ization constant co,o- Notice that 



(*o|*o) = |co,o|'(0|e "/"'^"-^^e "/"^^"-'^|0) 



APPENDIX C: UNIVERSAL SUPERFLUID 
PHYSICS 

In the present appendix we review and emphasize some 
of the hydrodynamic and universal properties of the su- 
perfluid phase as realized in Helium and BECs. 



(vI/o|vI'o) = |co,o|'E(-l)"^'"( — )" 



7 'm' 



and furthermore 



a^_-|0) = x/T|l,0> {al.f\Q) ^ VTx/2|2,0> 



n\\Q,n) 



1. Rigidity and Superfluidity 

We first review the Landau argument [ll| illustrating 
how the lack of rigidity causes dissipation. Consider a 
Bose fluid at zero temperature (a Bose condensate) with 
excitation spectrum given by ep. Suppose the fluid is 
moving with velocity v and interacts with something ex- 
ternal (what does it interact with is irrelevant, but one 
could think that it is a surface over which it is moving). 

The interaction can exchange energy and momentum 
with the fluid. The fluid has initial energy Mv^ /2 with 
M the total mass of the fluid. The interaction creates 
a quasiparticle with energy ep in the frame of the Bose 
fluid. In the lab frame the particle has energy 



ep= ep + p-v 



(CI) 



,,,|0)=0 a_j:\n) = ^\n-l) 



{a_£f\n) = V^V^T^ln - 2) (a_£)"|0) = 
Using the expressions above we find 



\a2a-~{a^^-{al.r\0) = n!(0|a^a:-|n, n) 

.M2 



|0) 



(n!)2(0|0) = (n!) 



n2 



which finally yields 



(M'ol^o) = |co,ol'E rdr Sn.min'-f 



n!m! 



by Galilean invariance since the velocity of the particle 
in the fiuid frame is vp = dep/dp and it \s vp + v in the 
lab frame. 

Energy will be exchanged if the fluid can lower its en- 
ergy. The fluid final energy is 



+ ep + p- V 



(C2) 



The condition for dissipation is thus ep + p- v < 0. The 
easiest way to satisfy the condition is when p is antipar- 
allel to V. We then obtain dissipation for 



P 



(C3) 



We see that for all dispersions softer than sound the min- 
imum on the right hand side of the last equation is zero 
and there is dissipation for any velocity of the fluid. If 
the dispersion is at least as stiff as sound, for velocities 
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smaller than some critical velocity there is no dissipa- 
tion and hence superflow. Since as long as the system 
is a quantmn Bose fluid, i.e. there exists a condensate, 
repulsion produces sound excitations as the elementary 
low energy quasiparticles of the system, the interactions 
provide the necessary rigidity for superfluidity. The free 
Bose gas, despite the macroscopic coherent occupation of 
the lowest momentum state, is not a superfluid as it is 
arbitrarily easy to destroy its coherence no matter how 
weak the perturbation. In fact, at energy scales in which 
the superfluid behaves as a free Bose gas, it dissipates 
and it is not super at all. This happens for wavelengths 
shorter than the coherence length of the superfluid as 
explained in a subsequent section. 



2. Scattering 

The purpose of this section is to illustrate the prin- 
ciples of superfluidity and of macroscopic exactness by 
an example. Among the experimental efforts on BECs 
figure studies of spontaneous emission and light scat- 
tering of these superfluid systems by W. Ketterle and 
collaborators l2^. They found that light scattering off 
BECs is suppressed to negligible values at long wave- 
lengths. This seems surprising at first since one can 
naively expect the typical n + 1 bosonic enhancement 
as observed in their emission experiments [26]. They have 
attributed the observed suppression of light scatte ring 
in BECs to destructive quantum interference effect s[l4l|. 
While their argument is certainly true, it is the purpose 
of this section to point out that the more universal reason 
that makes their argument, and hence their results, true 
is the Principle of Superfluidity. We show that the sup- 
pression of light scattering off a Bose Einstein Conden- 
sate is equivalent to the Landau argument for superflu- 
idity and thus is a consequence of the Principle of Super- 
fluidity. This means that the superfluidity is ultimately 
the reason for suppressed scattering at low wavelengths. 
We have previously presented these arguments in a short 
notep^. 

We will study next scattering processes in the super- 
fluid. In a scattering event, an incoming particle couples 
to the density of the system. The effective scattering 
Hamiltonian is 



d^rV{r, t)p(r) + H.C. 



where by Fourier transforming 



(C4) 



(C5) 



Hence 



ni^J^ ^9 / e-*(« '-'^*V(r) + H.C. (C6) 

is the interaction matrix element between the final 
and initial states of the external particle. Now, pg- = 



/ d^r p{r)e and at an arbitrary point i, the density 
is given by p{ri) = a|aj, where 



Oi = > e '''''"a^ a. 



thus 



= E^ 



k' ^ 



(C8) 



k,k' 



d^r p(r)e-'« '' 



k,k' 



k,k' 



and 



(CIO) 



where the energy and momentum transfered in the scat- 
tering event are 



Qf — qi and lo — uji ~ ujf 



(Cll) 



and where we have concentrated only on one q. Since at 
first order the amplitudes for the different momenta add, 
it is enough to do the calculations for one q. This will 
not be true for higher orders. Notice also that in the case 
of elastic scattering u; = 0. 

We need to insure adiabaticity in the interaction. Oth- 
erwise, a sudden hit on the system caused by an abrupt 
appearance of the scattering particle will do violent 
things to the system and cause it to generate entropy and 
heat up. We achieve this by the introduction of an extra 
term in the exponential of the interaction Hamiltonian[^. 
The interaction will be turned on at i = — oo and left on 
until a time t. At the end of the calculations, we will 
take the limit 77 ^ 



k+q k 



H.C. 



(C12) 



Calculating the time evolution of the Bose system in 
the presence of the scattering particle by means of time 
dependent perturbation theory we obtain for the ampli- 
tude 



J — 00 



n - 

(C13) 

where in our case the initial state is the ground state: 
l^i) = 1^0), Ei = 0. For the scattering amplitude we are 
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interested in to 7^ i, where 

J — 00 



fc+<? 



J(Em-En+huj-ihri)t/h 



(C14) 



and thus 



Em — huj — ihr] 



(C15) 



The scattering transition amphtude is, after separating 
the fc = 0, —q momentum states to account for the con- 
densate 



Cm{t) 



l/^\/Mo((*™|at|4'o) + (*m|a_d*o» 



Em + hoj — ihr] 



_^i{Ej-rL—hLO — ihri)t/h 



Em — hcu — ihr] 



Em + huj ~ ihr] 



^i{Ern+huj —ih'r])t / h 



Let us look more closely at the first two terms. They 
represent processes in which the external particle hits 
the atoms in the ground state and tries to excite them, 
i.e. knock an atom into or knock it out of the conden- 
sate. On the other hand, the entanglement in the ground 
state makes sound the only possible excitation at low 
momenta. In light of this, it is not too surprising that 
scattering is suppressed. It is just the principle of super- 
fluidity at work. 

Since the scattering particles couple to the density, the 
amount of scattering is proportional to the density com- 
ponent with momentum q, the so called density response 
function 



E 



Ejn — huj — ihr] 



(CI6) 



Notice from ((T7|) that 



at = u^bi - Vq'b-q' = u^bq - v^b^_^ (C18) 



Using these relations together with (|C8|) and (|C9|I we see 
that 



(*oiPd*™> = E(*oi4+/£i*' 



(C19) 



(Pq) = E Cm(*0|pd*m> + E <n(*m|pd*0> (C17) whcrC 



k+q 



(*olK^+,-&U,-- Vj:+^^b_j:_^){uf:bf: - vj:bl.)\^m) = ~Vf:^-Uf:{^o\b_j:_^J>j:\^m) 

= -^£+f"fc(*o|fo_fc_,A^|*™> + Vf:^.v^,{^omq) + bl-b_j:_^\^m) 
- v~,^.vj:{^o\S{qWm) = -vUqU^{^o\b_]:_J>]:\^m) (C20) 
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Since we are considering scattering with nonzero In this case the only final state not killing the matrix 

5{q) — Q. The case where q = corresponds to no mo- element is 



mentum transfer, i.e. no scattering and we discard it. 
Then the only final state giving a nonzero matrix ele- 
ment is 



thus giving 



(C21) 



(C22) 



For fc = we have similarly 

= -V^Wg-(«'o|&-g|*m) 



(C23) 



which gives 

The same calculation for k = —q yields 

(^olPfl*™) = VM^(*o|a-g|*m) 

= VMoUg-(*o|&-g|*m) 



(C24) 



(C25) 



(C26) 



with the same final state (jC24|) . Continuing the same 
type of calculation we also find 



J 



(*™|aL„-afcl*o> 



k+q 



I 



-9 -fe' 



^?^„-«?(*™|6t^/_,,|*o> (C27) 



Again, since g 7^ in our considerations, 5{q) — and 
we find as the only possible final state 



corresponding to a matrix element 



For fc = 



(*ml4l*0> = (*mlK-*l--l'9-6-g-)|*o) 



with final state 



and matrix element 



*m|4l*0) = Uff 



(C28) 



(C29) 



(C30) 
(C31) 
(C32) 



-^l-u; -^q 

Finally, setting fc = —q we obtain 

= -t;,-(^'„|6p«'o) (C33) 
whose final state is the same state ljC31|l 



Substituting l(n2T|l and (|n24ll into (|U2)l we find the 
two scattering amplitudes Cma (which corresponds to 
both fc = and k = —q since both have the same 
final state) and Cmi (corresponding to processes with 
fc 7^ 0, —q) to be 



^ ~Kf\/M^(Mg - Vq) ^i{Eff-hLj-ih-n)t/h 
Eq' — huj — ihrj 



(C34) 



y k k+q 



(C35) 

On the other hand, substituting ljC28|l and ljC31|) into 
(|C 2|l gives for the scattering amplitudes Cm' and Cm' 



j-VMo(ug-- Vq) ^i{Ef;+hu-ihn)t/h 
+ ftu; — ihrj 



(C36) 



Vffvrur , - 

q k k+q 



_^i{Ei:^.+ E^+huj-ihv)t/h 



hence 



-VqVM^{Uq- Vq) i{Es+huj+ihri)t / h 



E^ + hio + ihri 



(C37) 



(C38) 
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Plugging all these results into ljC17|) we obtain 



'J k k+q 



(C39) 



VqMo{u^ v^Y i(Eff-huj-ihn)t/h -iEfft/h _ ^gMpju^ v^Y -i(Eff+huj+ihr,)t/h iEfft/h 



E 



ihri 

E E 



k^O,-ql^O,-q 



^^k + q+^k-^-T'^'l^ 



+ So; + ifvq 

AEug+E^-hu,-thn)t/h -t{E^^^+E^)t/h 



E 



E E 



(C40) 



The higher momenta terms (fc 7^ 0, — in the density 
response function represent processes that are forbidden 
by kinematics constraints, hence they are dropped. In 
the case of light scattering conservation of energy gives 



cqi = cqf + 2csk 
while conservation of momentum gives 
Qt ^ qf + k- k 



(C41) 



(C42) 



thus 



Eq 



oUo (1 



mCs 
hq 



n-'q-' 



_ n?q'^ 

2^^ 

8m2c2 



2„2 



h^q 



(C47) 



for there is a particle excited with momentum k and an- 
other one with momentum —k. We have 



2 \E^ 



Eq 



- 1 



(C48) 



cqi 



If 
cqf 



(C43) 



which can only be satisfied simultaneously with the con- 
servation of energy condition in IIC42II by fc = 0. Thus 
the processes are forbidden. 

As for the multiplicative factor (ug*— v^)'^, it is an in- 
terference factor which provides suppression of light scat- 
tering as required by the necessary rigidity of the ground 
state. Expanding it out and using equations pt)|) . (|2(J|I . 
(ED), and (I23l we find, as g ^ 



(C44) 



Mo;7o(l + 



h^q^ 



2mMoUo' 



MoUoil 



2„2 



h^q 
2m?c; 



:) (C45) 



Hence 



2,2 

^q 



mCg 
hq 



1 + 



(C49) 




(C50) 



As can be clearly seen from HC30|) and (jC33p , repre- 
sents the probability of knocking a boson with momen- 
tum q out of the condensate, u| represents the probabil- 
ity of knocking into the condensate one of the boson that 
are correlated with momentum — (f, and is the inter- 
ference between the two processes. Putting everything 
together 



Eq 



hq 



MoUo 



h^q^ 



1/2 



hqcs 1 



h^q^ 
(C46) 



mCs 

hq 

hq 
2mCs 



1 



3h^q^ 



mCs 
hq 



(C51) 
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This goes to zero linearly as q ^ 0. The processes of 
knocking a particle with momentum q out of the corre- 
lated ground state and knocking a particle with momen- 
tum — g into the ground state both lead to the same final 
state. This is so because neither of these is an elementary 
excitation of the system, but a quantum superposition 
of them makes a Bogolyubov quasiparticle, which is an 
eigenstate of the system. Different processes leading to 
the same final state cannot be differentiated quantum 
mechanically and thus will interfere. The interference is 
destructive at small wavevectors and scattering is sup- 
pressed. Afore fundamentally, this result follows from 
the rigidity concomitant to the superfiuid state, which 
is ultimately due to repulsive interactions between the 
bosons. There is no scattering because the "Principle of 
Superfluidity" forbids it. These are universal properties 
valid for any superfiuid state which will hold exactly at 
low enough energy scales. 



APPENDIX D: SOME USEFUL RESULTS 

We present the calculation of some operators in the su- 
perfiuid ground state 
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APPENDIX E: EXPECTATION VALUE FOR 
NUMBER OF PARTICLES PER SITE IN THE 
MOTT GROUND STATE 



We conclude by emphasizing that if there is no repul- 
sion between the atoms, low energy excitations will not 
have a sound spectrum, t;^ = 0, = 1, leading to no su- 
perfiuidity and no suppressed light scattering. We thus 
see that suppressed light scattering follows from the rigid- 
ity concomitant to the superfiuid state, characterized by 
vj: ^ 0, which is equivalent to a nonzero sound speed. 
This rigidity is ultimately due to repulsive interactions 
between the bosons. There is no scattering because the 
"Principle of Superfiuidity" forbids it. These are univer- 
sal properties valid for any superfiuid state which will 
hold exactly at low enough energy scales. 



The expectation value for the number of particles per 
lattice site in a Mott insulating phase can be easily found 
y noticing that 
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and taking advantage of the commutation relations 
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we finally find 
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APPENDIX F: IMPOSSIBILITY OF 
SPONTANEOUS SYMMETRY BREAKING IN 
LOWER DIMENSIONAL SYSTEMS 

In the present appendix we review how a continuous 



\ \/NI \/N\ j symmetry cannot be broken even at zero temperature in 

I one dimension. This is the quantum version of the so 

[O'li (.Qq) j called Mermin- Wagner theorem "2^ on melting of order 

^/]V! in one and two dimensions by thermal fluctuations. An 

yfj^ ^ 1 ^ earlier reference with this result is Phil Anderson's lec- 

ture notes on solids 5^ where the melting of order by 
quantum fluctuation is treated too. We consider the ef- 
fective low energy Hamiltonian for a macroscopic system 
with a broken symmetry ground state: 



(a,|a,) = (n,) = 1 (E6) 



where is the frequency, or energy, of the excitation, 
and 6t are the annihilation and creation operators of 
elementary excitations from the ground state. Since we 
are supposing that this is the effective Hamiltonian of a 
continuous broken symmetry ground state, Goldstone's 
theorem implies that ^ as fc ^ [sDl- We, of 
course, have in mind the Bogolyubov Hamiltonian, but 
this applies to any Hamiltonian having a continuous bro- 
On the other hand ken symmetry lowest energy state. 

This Hamiltonian can be written in terms of coordi- 
nates and momenta. The coordinates would represent 
{n1) — the density at wave vectors k of the broken symmetry 
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breaking. This is done through the definition 

Q~^ = -kz^h + b]:) (F2) 

Pf^ = '\l^^^l-h) (F3) 

which are, of course, canonically conjugate variables, that 
is, [Q^,Pg] = i{Q£) = 0- The Hamiltonian then is easily 
seen to be a collection of harmonic oscillators for each k: 




In the ground state, the oscillators are centered on av- 
erage {Q-^ = 0. On the other hand, the uncertainty 
principle provides for density fluctuations ^ 0. The 
magnitude of the fluctuations can be estimated readily 
because the ground state energy is virialized. 



We thus see that ((5|) = 1/(2^^). Due to Goldstone's 
theorem, this quantity diverges at long wavelengths. This 
need not be a problem if the density of states vanishes 
at long wavelengths faster than the divergence. Other- 
wise, uncontrolled quantum mechanical density fluctua- 
tions will melt the ordered ground state and the symme- 
try will not be broken. This would lead to absence of 
superfluidity or phase stiffness of the ground state. In 
1-D, the total fluctuations is given by 



= j dk {QD = I ^ (F6) 



which for the Bogolyubov ground state, tOk k, diverges 
logarithmically at long wavelengths. Therefore, superflu- 
idity is not possible in one dimension. In the case that 
the finite size of the system cuts off the divergence, the 
system will still not be a superfluid, as this will happen 
when the energy spacings are experimentally discernible, 
which in turn means the system is not a even a fluid 
anymore for it does not behave like a continuum. 
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